Changes in tree functional composition across topographic gradients and through time in a tropical montane forest

Understanding variation in tree functional traits along topographic gradients and through time provides insights into the processes that will shape community composition and determine ecosystem functioning. In montane environments, complex topography is known to affect forest structure and composition, yet its role in determining trait composition, indices on community climatic tolerances, and responses to changing environmental conditions has not been fully explored. This study investigates how functional trait composition (characterized as community-weighted moments) and community climatic indices vary for the tree community as a whole and for its separate demographic components (i.e., dying, surviving, recruiting trees) over eight years in a topographically complex tropical Andean forest in southern Ecuador. We identified a strong influence of topography on functional composition and on species’ climatic optima, such that communities at lower topographic positions were dominated by acquisitive species adapted to both warmer and wetter conditions compared to communities at upper topographic positions which were dominated by conservative cold adapted species, possibly due to differences in soil conditions and hydrology. Forest functional and climatic composition remained stable through time; and we found limited evidence for trait-based responses to environmental change among demographic groups. Our findings confirm that fine-scale environmental conditions are a critical factor structuring plant communities in tropical forests, and suggest that slow environmental warming and community-based processes may promote short-term community functional stability. This study highlights the need to explore how diverse aspects of community trait composition vary in tropical montane forests, and to further investigate thresholds of forest response to environmental change.

Understanding variation in tree functional traits along topographic gradients and through time provides insights into the processes that will shape community composition and determine ecosystem functioning. In montane environments, complex topography is known to affect forest structure and composition, yet its role in determining trait composition, indices on community climatic tolerances, and responses to changing environmental conditions has not been fully explored. This study investigates how functional trait composition (characterized as community-weighted moments) and community climatic indices vary for the tree community as a whole and for its separate demographic components (i.e., dying, surviving, recruiting trees) over eight years in a topographically complex tropical Andean forest in southern Ecuador. We identified a strong influence of topography on functional composition and on species' climatic optima, such that communities at lower topographic positions were dominated by acquisitive species adapted to both warmer and wetter conditions compared to communities at upper topographic positions which were dominated by conservative cold adapted species, possibly due to differences in soil conditions and hydrology. Forest functional and climatic composition remained stable through time; and we found limited evidence for trait-based responses to environmental change among demographic groups. Our findings confirm that fine-scale environmental conditions are a critical factor structuring plant communities in tropical forests, and suggest that slow environmental warming and community-based processes may promote short-term community functional stability. This study highlights the need to explore how diverse aspects of community trait composition vary in tropical montane forests, and to further investigate thresholds of forest response to environmental change.

Introduction
Climate change is expected to drive major shifts in tropical plant communities [1][2][3][4]. Given the tight relationship between plant traits and environmental variation [5,6], understanding trait-based responses to climate may provide valuable insights into the processes that will shape the composition and functioning of plant communities under changing environmental conditions [7,8]. Tropical montane ecosystems are critical repositories of biodiversity [9, 10] and deliver many valuable ecosystem services [11]. Unfortunately, tropical montane systems are also believed to be especially sensitive to ongoing environmental alterations [12]. Changing geographic distributions is one commonly observed response of montane species to warming [13,14]. Indeed, many Neotropical montane forests are already exhibiting increasing relative abundances of tree species that were formerly distributed at lower elevations, presumably because these lower-elevation species are better adapted to hotter temperatures (i.e., are more thermophilic) than their high-elevation counterparts, and thus are able to increase in abundance or shift their ranges upslope (via increased recruitment at the leading edge and/or increased mortality at the trailing one) as temperatures increase [15][16][17][18].
In tropical montane forests, the functional composition of trees varies along elevation gradients. In general, tree species characterized by more conservative strategies for resource acquisition and use (e.g., higher leaf toughness, higher wood density [WSG], lower specific leaf area [SLA], and lower leaf nutrient concentrations) are more abundant at higher elevations, and species with more acquisitive strategies (e.g., higher leaf nutrient concentrations, higher SLA, and lower WSG) prevail at lower elevations [5]. These patterns are typically attributed to the steady adiabatic decrease of temperature with elevation [5] along with decreasing nutrient availability [5,19]. As such, anthropogenic global warming is expected to push the functional composition of tropical mountain plant communities towards greater relative abundances of species with more acquisitive traits [20]. In other words, changing species distributions due to warming have the potential to modify not only the floristic composition of communities, but also their functional composition (i.e., the relative abundances of acquisitive vs. conservative species) [20,21], which in turn can affect community dynamics and ecosystem functioning.
Changing precipitation rates can also have profound effects on plant communities [22,23]. Along tropical elevation gradients, the influence of precipitation on tree functional composition appears to be weaker than that of temperature [5]. This may be due in part to the fact that precipitation rates do not change consistently or predictably across elevations [24], as well as to the difficulty of translating precipitation to biologically meaningful metrics of plant water availability. That said, since climate change is affecting temperature and precipitation in these ecosystems [25,26], both factors clearly need to be considered if we hope to understand the responses of tropical tree species and communities to climate change.
While the abundance of low-elevation, thermophilic tree species is increasing in most Andean forests [16][17][18], the rates of change are highly variable [18]. This suggests that some biotic and abiotic factors prevent range shifts in certain plant communities [18,20,27].
Topographic variation is a critical driver of structure, composition and dynamics in tropical forests [28][29][30][31]. At a given elevation, topography shapes plant communities through a variety of mechanisms, including through effects on hydrology, nutrient availability, temperature, wind exposure, and biological interactions [32]. Bottomland and lower slope areas usually have greater soil moisture and nutrient concentrations than the upper slopes and ridges [33][34][35]. In addition, valleys and lower topographic positions often experience cooler temperatures relative to upper positions due to cold air drainage and reduced sun exposure [36,37], although the opposite patterns [38,39] or no trends of temperature variation from lower to upper topographic positions have also been reported for some tropical montane forests [35].
These differences in resource availability and micro-environmental conditions can have major effects on taxonomic and functional trait composition such that tree communities in lower  topographic positions typically include more species with acquisitive traits, whereas upper  positions and ridges are usually dominated by species with conservative traits [29, 34, 40]. However, it remains unknown how these patterns of plant trait distribution driven by local topographic variation relate to the species climatic preferences considering regional gradients of elevation.
At a particular site, topographic variation may modulate forest responses to climate change through various mechanisms. Topography may strongly regulate small-scale climate variation and plant community composition in montane environments [38,41,42], prevent some species from tracking climates [18,27], and define the location of climatic refugia for some species [43][44][45]. Under increasing temperatures, the tree communities in lower topographic positions that harbor acquisitive species with a greater potential for fast resource uptake and investment [46,47] may be expected to respond more-strongly in terms of growth than communities at upper positions that have more conservative species. In contrast, higher temperatures could promote the colonization of more acquisitive species into habitats at upper positions, and thereby affect the floristic composition of these communities. Thus, trait-based responses of different demographic groups can indicate the mechanisms that determine plant community responses to climate change along topographic gradients. Particularly marked differences in trait composition may occur between mortality and recruitment demographic groups, as dying trees were established under environmental conditions likely different than present ones, under which new recruits are succeeding.
We investigated how functional trait composition and community climatic indices vary across spatial and temporal gradients in a topographically complex tropical montane forest in the Andes of Ecuador. Under the premises put forward by Trait Driver Theory that sustained trait-based responses to directional environmental variation can be described by communityweighted means (CWM) [48], we focused our analyses on exploring how CWM vary in response to environmental gradients, through time and across tree demographic groups. Following the same theoretical context, to have a broader understanding of trait distribution within the studied tree communities, we explored how three additional community-weighted moments: variance (CWV), kurtosis (CWK) and skewness (CWS) varied across topographic gradients and time [48]. Using data collected over an eight-year period, we tested the hypotheses that 1) topography affects tree community trait composition from lower to upper topographic positions, which may reflect on community climatic indices and the relative abundance of tree species with acquisitive vs. conservative traits, and 2) community thermal and precipitation preferences may be increasing through time, such that communities will exhibit increasing relative abundance of species with acquisitive strategies and with higher thermal optima, due to the upslope shifts in the ranges of species. Understanding the role of topography on tree functional trait distribution and climatic optima is critically needed to help explain the observed variation of tropical montane forest sensitivity to climate change.

Results
In total, we measured 1594 dicotyledonous trees in our 18 plots-1341 individuals in the first census, and 1315 individuals in the last census, 1062 trees stayed alive over the eight-year census period, 279 individuals died, and 253 individuals were recruited. The average annual turnover rate for individuals was 2.87 (± 0.29) %. Of the measured trees, 1425 (89.4%) were identified to species, 105 (6.6%) were identified to genus, 49 (3.0%) were identified to family level. Fifteen individuals (1.0%) remained undetermined. In total, we identified 220 species or morphospecies.

Trait composition, community climate indices, and topography
Topography significantly affected the community trait composition (CWMs) of eight functional traits ( Table 1, Figs 2 and 3, S1 and S2 Tables). At lower TPI values, the tree communities had high average bark thickness, leaf N and P content, leaf area (LA), specific leaf area (SLA), sapwood-specific conductivity (KS), vessel diameter, and low average values of leaf toughness, vessel density and wood density (WSG). At higher TPI values (towards upper positions), this pattern of trait composition reversed, reflecting the expected gradient of acquisitive to conservative strategies of plants along the topographic gradient.
In our data set, the functional traits and climatic optima of species were correlated. Leaf N and P concentrations, SLA, and vessel diameter correlated positively with species' optimal mean annual temperature, whereas leaf toughness and vessel density correlated negatively with species' optimal mean annual temperature (  correlated in a positive and negative fashion, respectively, with species' mean total annual precipitation ( Table 2).
The repeated linear mixed effect models indicated that CWVs, CWKs and CWSs were less affected by topographic variation than CWMs (S4 Table and

Fig 2. Scaled Community Weighted Means (CWMs) of leaf and stem functional traits as a function of topographic variation Topographic Position
Index (TPI). All traits except for bark thickness and sapwood-specific conductivity varied significantly with TPI (S2 Table). The slopes of the linear regression lines represent the first (2008) and last (2016) sampling years were not significantly different for any trait. Each point represents a forest monitoring plot.
https://doi.org/10.1371/journal.pone.0263508.g002 Community-weighted skewness (CWS) decreased for leaf P concentration toward areas with higher TPI values. Community-weighted kurtosis (CWK) was not significantly affected by topography. For our study plots Community Thermal Index (CTI) significantly decreased with TPI. This indicated that plots at lower topographic positions had higher relative abundance of species with hotter thermal optima compared to the upper slope forests (Linear regressions: CTI yr1 = 18.91-1.52 x TPI, R 2 = 0.57, P< 0.001, CTI yr8 = 18.67-1.36 x TPI, R 2 = 0.54, P <0.001; Figs 3 and 4A, S2 Table). In fact, plots located towards lower topographic positions had CTI's up to 4˚C higher than those at upper topographic positions. Similarly, Community Precipitation Index (CPI) decreased with TPI, such that plots at lower topographic positions had greater relative abundances of species adapted to wetter environments (higher CPI) compared to plots located at upper topographic positions (Linear regressions: CPI yr1 = 1987.65-120.46 x TPI, Table 2. Non-parametric correlations between species functional traits, and species climatic optima, and number of species included in the respective analyses. Statistically significant values are presented in bold.  Table). Plots at lower topographic positions had CPI values that were up to 400 mm greater than those of communities located at upper topographic positions. Community-weighted skewness (CWS) increased along topographic position for CTI and CPI (S4 Table and S2 Fig).

Community functional trait composition and climate indices through time
The CWMs did not change over the study period along topographic gradients for any of the ten functional traits that we measured (Fig 2, S2 Table). Similarly, CWMs of the demographic components of the tree community (i.e., dying, surviving and recruiting trees) had statistically significant different mean values in only one functional trait: vessel diameter (S5 Table and S3  Fig). Specifically, CWMs of stem vessel diameter of trees that stayed alive through the study period was significantly higher than of the community of recruits (Tukey HSD P<0.05).
Our study plots did not experience directional changes of their Community Temperature Index (CTI) over the eight-year study period (TRplot mean = -0.02 ± 0.01 per year, binomial probability distribution P = 0.759). The mean CTI of the entire plant community did not change over time (Fig 4A, S2 Table), nor across the three demographic groups (Fig 5A), as they were only affected by TPI following the same trend of the entire plant community (S5 Table).
As with CTI, the Community Precipitation Index (CPI) of our plots did not experience directional shifts over time (PRplot mean = -1.26 ± 1.26 per year, binomial probability distribution P = 0.240), and their mean CPI did not change over time (Fig 4B, S2 Table). Demographic groups did not differ in their CPI either but responded to TPI as observed for the entire plant community (Fig 5B, S5 Table and S3 Fig). None of the three community-weighted moments (CWK, CWS, CWV) changed significantly over time in our forest monitoring plots (S4 Table and

Discussion
Our results indicated that the topographic position of our plots strongly influenced the floristic and functional composition of the tree community and revealed unrecognized links between topography and species climatic optima in tropical montane forests. Furthermore, our censuses and demographic analyses indicated that forest functional and climatic composition remained relatively stable at our study site over the eight-year study period. As discussed below, we propose that variation in soil conditions rather than temperature mediated the observed patterns of community trait composition and that functional stability may be driven by various factors, including high species diversity, the location of our study site at an ecotone, low rates of climatic warming, and topographic complexity.
In agreement with our first hypothesis, the observed gradient in CWMs of acquisitive to conservative traits in plots at lower-to upper-topographic positions was in line with previous research in tropical lowland and montane forests (Figs 2 and 3, S2 Table) [28-31, 40, 49] and supports the pattern of fine root traits found in the same study plots [50]. The lower values of CWVs toward high topographic positions for four of the examined traits were also in line with the trends observed for CWMs, and confirm the role of topography as an environmental filter [48]. It is interesting to note that functional trait composition measured as CWMs along topographic gradients (from lower-to upper positions) mirrored the patterns observed along large elevation gradients in tropical forests [5], emphasizing the role that local factors, and not just large-scale climate gradients, can play in structuring trait distribution in these ecosystems.
An unanticipated finding of this study was the strong relation between topographic variation and community climatic indices. Indeed, CTIs in plots at lower topographic positions were up to 4˚C higher than in plots at upper topographic positions, which would correspond to a difference of~700 m in elevation assuming a lapse rate of 5.5˚C per 1000 m (in contrast to the 1˚C range of CTI expected across our study plots which are separated by a maximum of just 176 m, Figs 3 and 4A). The higher CTI of forests in lower topographic positions contrasts with the idea that lower topographic positions and valleys experience decreased temperatures due to lower solar energy inputs, cold-air drainage and high canopy cover [36,37]. For tropical humid montane forests, the only study we know of that reports field-recorded temperatures along topographic gradients did not find statistically significant differences between lower and upper topographic positions [35], suggesting that other factors (e.g., lower variation in aspect, wind) may homogenize environmental temperatures in these ecosystems. In consequence, although topography may have affected microclimatic conditions at our site, we think it is unlikely that these climate differences could have caused the gradients of up to 4˚C in CTI and 400 mm in CPI, especially given the small differences in elevations between sites (Figs 3 and  4B). Rather, we hypothesize that topography is a filter for plant functional traits that correlate with the species' thermal and precipitation optima. In particular, the effects of topography on local soil conditions and hydrology may influence the establishment of plant communities with varying requirements of water availability [22], which in turn relate to species' thermal preferences. Indeed, seven of the ten functional traits analyzed had statistically significant correlations with species optimal temperature, and two with optimal precipitation (Tables 1 and  2).
The effects of topography on climatic composition may have been particularly strong at our study site, which was located at the upper elevational limit of lower montane forests [51]. In this transition zone, topographic complexity could facilitate the establishment of lower montane species at their maximum elevations, especially in lower topographic positions that are relatively sheltered and that have higher soil moisture and nutrient content. In contrast, upper topographic positions that are more exposed with poorer and periodically drier soils may be inhospitable for the lower-elevation species and favor stress-tolerant species from higher elevations that are less thermophilic and have lower precipitation optima. In fact, this process contributes to understand the positive effect of TPI on CWSs in the two community climatic indices (S4 Table and S2 Fig). For both of these (CTI and CPI), skewness became more positive (i.e., higher densities of low CTI and CPI values, left leaning distribution) toward high topographic positions, indicating that as TPI increased, the representation of species with higher thermal and precipitation optima decreased. Thereby, these results highlight the possibility that in montane forests topography may be a particularly powerful factor for structuring plant communities around ecotones, magnifying stress gradients and creating mosaics of species with different thermal and precipitation optima at fixed elevations, thereby enhancing community complexity at landscape scales.
Stable trait composition in spite of the observed floristic change (Fig 1B) probably occurred because the recruiting species were functionally equivalent to the species that died, suggesting that environmental pressures were not acting selectively for or against species with certain functional traits. Rather, these results suggest that we witnessed the dynamics of a diverse mature forest where tree species turnover did not cause directional changes of functional composition, possibly due to high species diversity and functional redundancy [52][53][54]. Functional redundancy is expected to enhance the resilience of ecosystems under changing environmental conditions [55], and may be buffering directional changes of primary productivity in Neotropical montane forests where CTI is increasing due to rising temperatures [56]. The fact that our individuals turnover rates fall within the range of those estimated for humid tropical Andean forests [57] also suggests unaltered natural forest dynamics.
High rates of species turnover at upper and lower TPI sites were most likely not a result of higher environmental temperatures as we did not detect increasing CTIs (Fig 4A) or higher representation of species with acquisitive traits (Fig 2). Rather, the rapidly increasing rates of nutrient deposition in our study area [58,59] could have disadvantaged species adapted to poor soils and to a lesser extent facilitated the colonization of species adapted to more fertile soils [60]. At lower topographic positions, higher rates of floristic change may be a result of higher rates of turnover of individual trees compared to upper positions [29,30,58].
It is noteworthy, that the only trait difference between demographic groups that we found was the larger vessel diameter in surviving trees vs. recruited, which suggests that recruiting individuals are of species that sacrificed rapid resource uptake for protection from water deficits that may cause embolism [61]. Our study area experiences short dry spells every few years, which may affect tree recruitment; but it is unknown whether the frequency and intensity of these events is shifting due to climate change. More appropriate comparisons across demographic groups will require longer monitoring periods in which the recruiting trees have experienced a fuller range of biotic and abiotic constraints that ultimately determine the species composition of the entire plant community, or comparisons of the composition of the community demographic groups between forests censuses spanning shorter periods.
Finally, rates of climate warming in Andean forests ecosystems are predicted to be highly variable across elevations [26], potentially contributing to the heterogeneous rates of migration and functional change observed in forests throughout the Andes [18]. In this context, our study site seems to have experienced comparatively low levels of climate change and environmental forcing during our study period, perhaps below a response threshold, resulting in stable community climate indexes and functional trait composition especially given that our study was conducted over just an eight-year monitoring period. Indeed, the estimated annual temperature increase of 0.013ºC for our study area [62,63] lies well below the Andes-wide average warming rate of 0.6ºC [18]. Finally, complex topography has been demonstrated to contribute to lower rates of warming in montane areas [37,64]. Therefore, a combination of biotic and abiotic factors may contribute to functional stability in our study sites.
Micro-environmental conditions are increasingly recognized as critical drivers of plant responses to climate change [44, [65][66][67]. Indeed, in montane environments, fine-scale environmental heterogeneity may mediate the reorganization of plant communities under climate change in ways that do not necessitate shifts in species' altitudinal distributions, but rather follow specific combinations of topographic and micro-environmental conditions [41]. This notion emphasizes the need of intensive and long-term vegetation and micro-climate monitoring at landscape-scales to explore multidimensional forest responses that depart from upslope and poleward expectations of species movements under climate change [45,68], especially in topographically complex and species-rich tropical montane forests.
While our eight-year monitoring period is a relatively small window of time over which to document forest responses to climate change, previous research has found changes in the community climatic indices in Andean forests over periods as short as four years [17]. Overall, slow environmental warming, community-based processes, fine-scale environmental conditions, and the location of our plots within a transition between forest types all appear to have resulted in community functional stability. In a broader sense, our findings highlight the need to investigate forest response thresholds and how factors beyond climate affect tree trait distributions in tropical montane forests and present a cautionary call for studies exploring trait composition along elevation gradients as they may suffer from systematic topographic sampling biases at different altitudes.

Study site
This study was conducted in montane forests on the eastern slopes of the Andes of southern Ecuador (Reserva San Francisco, 3˚58'S, 79˚04'W) in a transition area between lower and upper evergreen montane forests [51]. The area has a mean annual temperature (MAT) of 15C and receives a mean total annual precipitation (MAP) of 2100-2200 mm. The area has been warming at the rate of 0.13ºC per decade since the early 1960s, but there have been no directional changes of precipitation [62,63]. The atmospheric deposition of nitrogen and phosphorus has increased during the last decades [58,59].
The study area has a complex topography with slopes of varying degrees of steepness. Soils are heterogeneous but are generally poor in nutrients [29,30,69,70]. Here, as in other tropical forests, topography strongly affects soil nutrient concentrations and hydrology, with soils at lower topographic positions exhibiting substantially higher concentrations of nitrogen, phosphorus and micronutrients, which affects species composition, tree trait composition, forest turnover and productivity rates [29,30,69,71].

Experimental design and field measurements
In 2008, eighteen permanent plots measuring 20 x 20 m each were established in a narrow altitudinal belt of < 180 m of elevation, between 1900 and 2100 m a.s.l. In the plots, we marked, identified, and recorded the diameter of all dicotyledonous trees �5 cm dbh (diameter at breast height). In 2016, we re-measured the trees in the plots, documenting tree death and recruitment.
The plots were distributed with a minimum separation distance of 50 m to capture topographic variation. For each plot we obtained a Topographic Position Index (TPI) derived from the relative position of the plot's corresponding raster cell in a digital elevation model (DEM) related to the average elevation of the surrounding cells [72]. We used a DEM with 10 m resolution generated from stereo aerial photos by aerotriangulation [73] to calculate the TPI for a circular neighborhood around each plot with a radius of 200 m [74].

Tree functional traits
For 158 of the 200 species present in our monitoring plots, we measured ten leaf and stem traits expected to respond to topographic variation as indicated in Table 1. Trait measurements were taken between 2008 and 2014 from trees distributed in our study area, using standard methods as described in S6 Table. For individuals identified only to their family or genus, we used mean trait values of their lower taxonomic levels.

Species climatic optima
We estimated the thermal and precipitation optima of the tree species present in our study plots based on the locations of observation and collection records relative to large-scale climate patterns following protocols modified from Fadrique et al. (2018) [18]. Specifically, we downloaded all available georeferenced records of the target species from the Andean countries of Venezuela, Colombia, Ecuador, Peru, Bolivia and Argentina from the Botanical Information and Ecological Network (BIEN) database [75]. We added distribution data from the BioWeb database [76] for those species that were absent or that had <20 records in BIEN. To minimize possible bias, we eliminated records with obvious georeferencing errors or that fell outside the Andean region and we only used one record per coordinate for each species to avoid including duplicates. Next, we extracted the mean annual temperature (MAT) and mean total annual precipitation (MAP) values at the record coordinates from the CHELSA extrapolated climate map [77] with a spatial resolution of 30 arcsec (approximately 1 km 2 at the equator). We then calculated each species' thermal optimum as the average of the extracted MAT values, and each species' precipitation optimum as the average of the extracted MAP values. We obtained the thermal and precipitation optima for 154 (70% of the species pool) species present in our plots, representing on average 70% of the plots basal area (S7 Table), using a mean of 187 observations per species.

Functional trait composition and community climatic indices
We used standardized trait values (mean = 0, standard deviation = 1) to calculate plot level community-weighted moments, a measure of community trait composition in which species' trait values were weighted by their basal area in a plot (i.e., basal area = summed cross-sectional stem area of all conspecifics in a plot). We calculated community-weighted moments for the first and last year of our monitoring period to evaluate directional changes in overall community trait composition following the equations A16, A17 and A18 described in Enquist et al. 2015 [48]. In this way, we captured variation in community-level trait composition (i.e., CWMs), but also on community-level trait distribution (i.e., CWVs, CWSs, CWKs).
Community-weighted means (CWMs) were also calculated for three demographic groups: trees that died between censuses, trees that survived, and trees that were recruited during the census period to investigate the associations between these demographic processes and functional trait composition [78].
As we did with species functional traits, we calculated the community-weighted moments of the species' thermal and precipitation optima in each plot (with species weighted by their plot basal area) to generate a community temperature index (CTI) and a community precipitation index (CPI) for each plot. We calculated the CTI and CPI for each plot for the first and last years of the monitoring period, and for the separate demographic subsets (i.e., dead, recruiting, surviving trees). For comparison, we also calculated CWMs using the number of stems per plot to measure trait representation in the community. Our results indicated that both ways of calculating CWMs yield highly correlated values (correlation coefficients between 0.73 and 0.99, and were highly statistically significant in most cases (S8 Table). We therefore conducted the statistical analyses using CWMs calculated with plot-scale basal area, because it better represents the contribution of the single species to aboveground biomass and ecosystem functioning.

Data analyses
In our dataset, elevation and TPI were significantly correlated (pairwise correlation = 0.72, P = 0.001). Thus, we performed various statistical analyses using PERMANOVAs and linear mixed models that included elevation and compared them using R 2 values and Akaike Information Criterion (AIC) [79] to evaluate the relative effect of elevation and topographic variation in our short topographic gradient. The approach and results of these analyses are presented in S9 and S10 Tables. Since these analyses consistently indicated that topographic variation was the main determinant of floristic and trait composition in our permanent plots, we performed subsequent statistical analyses excluding the elevation variable. We also tested for spatial autocorrelation in our trait measures to assess the need to include terms describing spatial distribution of the study plots in the statistical analyses [80]. There was no evidence suggesting spatial autocorrelation in our plot-scale trait measurements (see S11 Table for methods and results); we therefore developed further analyses without spatial correlation terms.
Floristic composition. To understand the effects of topography on species composition, we ran PERMANOVAs on Bray-Curtis species dissimilarity matrices using species plot scale basal area as a measure of species abundance, with 999 permutations in the 'vegan' statistical package [81] in R [82], and we used non-metric multidimensional scaling (NMDS) to visualize the results (S1 Fig We first explored the community-wide effects of TPI on species composition using data of our first monitoring year. Second, we evaluated floristic shifts as a function of time and topographic variation; thereby plot scale species composition was modeled as a function of time, TPI and their interaction. Finally, we tested for differences in species composition between dead and recruited trees, modeling species composition as a function of both demographic groups, TPI, and an interaction term.
We also calculated annual turnover rates for each plot which was estimated as the average of the annual percentage of recruitment, and annual percentage of mortality [83].
Community climate indices and trait composition. We used repeated measures linear mixed models to evaluate how trait composition measured as community-weighted moments (CWM, CWV, CWS, CWK) in each of our ten functional traits, and community climatic indices (CTI and CPI) changed across gradients of topography and through time. In our models, CWMs of a given functional trait or climatic index was a function of time, TPI, their interaction, and plot was included as a random factor. Similarly, we calculated repeated measures linear mixed models to evaluate the effect of topographic variation and time on the remaining three community-weighted moments (CWV, CWK, CWS). In these models, each community-weighted moment was modeled as a function of time and TPI, with plot included as a random factor. For statistical analyses we used the packages lme4 and lmerTest in R [84,85].
We also used differences between the initial and final censuses to obtain annualized rates of change in CTI and CPI per plot (TRplot and PRplot, respectively) [18]. We calculated binomial probability distributions with TRplot, and PRplot values against the null expectation of an equal number of plots decreasing and increasing in CTI/CPI over time.
Community climatic indices and trait composition across demographic groups. We used linear mixed models to assess differences in trait composition (only for communityweighted means CWMs) or climatic indices across demographic groups. Our models' fixed effects were demographic group (i.e., dead, surviving, and recruiting trees), TPI, and an interaction term; plot was included as a random factor. Then, Tukey-HSD was used to compare means among demographic groups when the interaction term was statistically significant.